Topological Entanglement in Abelian and non-Abelian Excitation Eigenstates 
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Entanglement in topological phases of matter has so far been investigated through the perspective 
of their ground-state wavefunctions. In contrast, we demonstrate that the excitations of Fractional 
Quantum Hall (FQH) systems also contain information to identify the system's topological order. 
Entanglement spectrum of the FQH quasihole (qh) excitations is shown to differentiate between the 
Conformal Field Theory (CFT) sectors, based on the relative position of the qh with respect to 
the entanglement cut. For Read-Rezayi model states, as well as Coulomb interaction eigenstates, 
the counting of the qh entanglement levels in the thermodynamic limit matches exactly the CFT 
counting, and sector changes occur as non-Abelian quasiholes successively cross the entanglement 
cut. 

PACS numbers: 63.22.-m, 87.10.-e,63.20.Pw 



Topologically ordered systems are not characterized by 
local order parameters; non-local concepts, such as quan- 
tum entanglement [1], have been extensively used in re- 
cent years to describe such phases of matter. The favorite 
method of analyzing the entanglement - entanglement 
entropy (or its topological part for gapped systems Q) - 
does not result in a unique characterization of the system: 
different states of matter can have identical entanglement 
entropy. Complicated topological phases, such as FQH 
states, are fully described by a multitude of universal pa- 
rameters, notably braiding matrices [3[, which are related 
to the properties of the FQH excitations under adiabatic 
exchanges in space- time. In finite systems, the braid- 
ing matrices are impossible to obtain and the question 
arises whether the universal properties of a topologically 
ordered state are obtainable via the entanglement of its 
excitations. Although scarcely addressed in the existing 
literature 0, [H[ , the question is pertinent also in view of 
the phases of matter that can only be distinguished by 
their excitation spectra Q. 

Recently, it was proposed [?J that the entanglement 
spectrum (ES), i.e. the (negative logarithm of the) full 
set of eigenvalues of the reduced density matrix p A is 
a rich source of information on the topological order in 
FQH ground states. Reduced density matrix p A of the 
subsystem A of a pure FQH state on the sphere 
(Fig. [1]) is given by the usual trace p A = Ttb\^)(^\ over 
the complementary subsystem B. The levels of p A can 
be classified according to the number of particles Na and 
orbitals Ia in A, as well as the z-axis projection of the 
angular momentum, Lf. The multiplicities and relative 
energy spectrum of p A matches that of the edge modes 
[7l-ll0j. For ground states of Coulomb Hamiltonians in 
the same universality class with a FQH model state, the 
ES typically displays a branch of low-lying (high prob- 
ability) levels, very similar to those of the model state, 
accompanied by spurious levels at high entanglement en- 



ergy (low probability). The gap between the low and high 
levels, properly defined by taking the "conformal limit" 
@, was conjectured and numerically substantiated to re- 
main constant upon increasing the system size. 




FIG. 1. FQH sphere with a monopole of magnetic flux 
quanta in the center, partitioned into hemispheres A and B, 
and containing the bosonic Moore-Read state with two qhs. 
We start with two separated qhs (left) and drag one qh from 
the north to the south hemisphere (middle). The moving qh 
is azimuthally delocalized. We end up with a qh twice the 
charge at the south pole (right). 

In this Letter we show that the ES of FQH excita- 
tion states contains information to identify the universal 
properties of topological phases of matter. We consider 
model wavefunctions, such as Laughlin Moore- Read 
[lH and Read-Rezayi [Hj], whose ground- and excited 
states with localized qhs can be expressed as Jack poly- 
nomials [3, [H|. Furthermore, we consider the eigen- 
states of Coulomb interaction potential with the impuri- 
ties that pin the qhs at a circle of latitude 

mini. The 

ES of a given FQH excitation is monitored as the qhs 
are moved across the cut (Fig. [Tj). This reveals that the 
ES of the excitations can probe different CFT fermion 
number sectors, that it gives the correct counting of the 
edge states in the thermodynamic limit, and that it is 
extremely sensitive (within a single magnetic length) to 
whether non-Abelian qhs are on the same or opposite 
sides of the entanglement cut. The latter property can 
be taken as a simple manifestation of the non-Abelian 
nature of the phase. These findings are corroborated in 
studies of the realistic Coulomb interaction eigenstates. 
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FIG. 2. (Color online) Conformal- limit ES of the Laughlin model state of N = 12 bosons with a single qh, localized at one of 
the poles (a), the equator (c) and in between (b). The location of the qh is given by the dot above each root partition. 



We also analyze the behavior of the quasielectron excita- 
tions, as well as the quasihole and quasielectron fluctua- 
tions when close to the entanglement cut. 
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FIG. 3. (Color online) Conformal-limit ES of the Moore- Read 
model state for N — 20 bosons with a unit flux added, (a) 
Abelian vortex 0202... 02. (b) Two fractionalized e/2 non- 
Abelian qhs. 
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FIG. 4. (Color online) Conformal-limit ES of the Z4 Read- 
Rezayi state of N — 28 bosons and a unit flux added. Abelian 
vortex (a) is fractionalized and the non- Abelian e/4 qhs are 
moved to the opposite pole one at a time ((b),(c)), causing 
the sector change. 

The orbital of an electron confined to the lowest Lan- 
dau level (LLL) moving on the Haldane sphere [HI 



(FigJTJ can be written as m (z) = A/" m £ m , where z — 
x+iy is the complex 2D electron coordinate, the quantum 
number m is the L z eigenvalue, and JV m is a geometri- 
cal normalization factor. "Conformal limit" [9[ is defined 
as the geometry where the normalization factors Af m are 
all equal to 1. General TV-electron LLL states are ana- 
lytic polynomials V F (zi, . . . , z N ) = H i<:j (zi -Zj)xV B , 
which can be factorized into a Vandermonde determi- 
nant and a bosonic wave function \£ b • We may therefore 
focus on the systems of charged bosons in the LLL. Ar- 
bitrary is expandable in terms of symmetric mono- 
mials indexed by a partition A represented by occupa- 
tion numbers n(A) = {n m (A),m = 0,1,..} of the or- 
bitals (j) m . Certain FQH wavefunctions are, however, 
defined by a single root partition - all the remaining 
partitions in the expansion of ^ b are derived from it 
via "squeezing" operations This includes the Read- 
Rezayi [13[ series of trial states for bosons at fill- 
ing factors v = N/N$ = k/2 which can be identified 
with a family of Jack polynomials (Jacks) J" , param- 
eterized by a = — (k + 1) and indexed by a partition 
rc(A fc ) = (fcOfcOfc...) i.e. ^ R (x ({*<}). Apart 



from FQH ground states, the Jacks also yield wavefunc- 
tions of the qh [l4| and quasielectron jl9[ excitations, 
created by varying the magnetic flux through the sphere. 
If the flux is increased by one unit, k non- Abelian qhs 
of charge e/k each are nucleated and pinned at loca- 
tions wi, . . . , Wk- Let us consider such an arrangement 
of the qhs where the first ri\ qhs are fixed at one pole 
of the sphere, ri2 = k — 1 — n\ are at the opposite pole 
and the remaining mobile qh is at Wk in between the 
poles. The wavefunction for states with a qh at Wk 
close to e.g. the north pole is given by a single Jack 
J ' 



-(fc+i) 



abab...ab 



({zi}) with the root partition abab . . . a6, where 
a = (k — 1 — An)/2, b = k — a and An = ri2 — ri\ [15|- As 
the qh at Wk is moved from the north to the south pole, 
the wavefunction mixes in several other Jacks 



14, Of 



and we can track its progression. When the qh reaches 
the south pole, the wavefunction is again a single Jack, 

.a+16-l d z i})- 



T-(fc + l) 

J a+lb-la+lb-1 
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(b) m=6 
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FIG. 5. (Color online) Conformal- limit ES of the Coulomb v — 1/2 state of N = 12 bosons with a single qh, localized using 
delta impurity potential, at one of the poles m — (a), the equator m — 12 (c) and in between m — 6 (b). The cut is defined 
by Na — 6 and l A 12. 



Jack qh wavefunction for k = 1 describes a single 
Abelian qh, localized at one of the orbitals of the sphere. 
In Fig. [2] we show the numerically calculated ES in the 
conformal limit for the Laughlin state of N = 12 bosons 
and a single localized qh. The cut is fixed such that the 
subsystem A contains 1^ = 12 orbitals and Na = 6 parti- 
cles. Above each plot in Figs. [2]|4j we give the root con- 
figuration used to generate the wavefunction and draw 
the cut in orbital space (vertical line), showing the high- 
probability ES levels in the inset. When the qh is at 
one of the poles (Fig. [2fa)), the level counting matches 
that of the U(l) chiral boson CFT. As we move the qh 
towards the equator (Fig. EJb)), the levels spread up- 
wards but the CFT counting remains unaltered. The 
spreading reflects the quantum-mechanical oscillation of 
the qh around the cut but is expected to be confined 
within a magnetic length around the equator; away from 
the equator, the counting of the qh entanglement lev- 
els remains a faithful representation of the FQH state. 
The CFT counting is lost only when the qh is situated 
exactly on the entanglement cut, which effectively splits 
the qh in two (Fig. Etc)). Using Jacks we obtained the 
ES of larger systems than those attainable by exact di- 
agonalization and N = 12 is presented only to facilitate 
comparison with the Coulomb case below. 

For k = 2 Moore- Read state, two non- Abelian qhs are 
formed when a unit flux is added. When both qhs are 
at the same pole, the wavefunction is represented by the 
root 0202 ... 02. When a single e/2 qh is fixed at one pole 
and the other is moved to the opposite pole, the Jack 
root is 1111... 11. ES for these two cases is shown in 
Fig. [3] and demonstrates that a sector change has taken 
place as the non- Abelian qh is moved across the cut: the 
counting 1,1,3... has changed into 1,2,4... The two 
different countings represent the two topological sectors 
of the theory, the even and odd fermion-number sectors 
respectively. Sector change occurs immediately as the 
qh crosses the entanglement cut. The thermodynamic- 
limit counting of the two excitation wavefunctions is in 
one to one correspondence to that of the excitations 
above the ground state for an even (20202000 . . .) and 



odd (20201000 . . .) number of particles. 

An interesting question is what happens when more 
than two topological sectors are present in the the- 
ory. The simplest example is the Z4 Read-Rezayi 
state where the three sectors' counting is given by 
the excitations above the ground states with N = 
0, 1, 2 (mod 4) number of particles: 4040404000 . . ., the 
01 sector 4040404010000 . . . (equivalent to the a[ sector 
40404040300000 . . .) and the a 2 sector 4040404020000 . . . 
(02 is its own conjugate). In Fig. 2] we show the ES 
for the Z4 state of N = 28 bosons with a single flux 
added. Starting from an Abelian vortex localized at one 
of the poles 0404 . . . 04 (Fig. [1(a)) with the counting 
1, 1, 3, 5, . . ., we transfer a single non- Abelian qh to the 
opposite pole. If we do this once (Fig. H(b)), we °b- 
tain the state with the root 1313 ... 13 and the count- 
ing 1,2,5,8...; doing it once again results in a root 
2222 ... 22 and the counting 1, 2, 6, 9, . . . (Fig. g^c)). 

Our main findings for the model states - the counting 
of the qh edge modes and the sector change upon crossing 
the cut - carry over to the realistic systems with LLL- 
projected Coulomb interaction. The problem of qh local- 
ization was solved generally in Ref. 16. For the Laughlin 
filling v = 1/2, we simulate the probe by a delta- function 
impurity potential [TtJ of weig ht 0.005/ V^/2, localized 
on the orbital m (Fig. [1]). This potential is weak enough 
not to cause any level crossing in the spectrum, but only 
splits the degenerate multiplets of states. In Fig. [5] we 
show the ES of the Coulomb ground state in the presence 
of a localized qh by superimposing it on the model state 
from the Fig. [2j A finite entanglement gap separates 
the CFT branch from the generic Coulomb continuum of 
higher entanglement-energy states. CFT branch, how- 
ever, displays the same counting as the model state. The 
Coulomb continuum largely remains fixed as we move the 
qh towards the equator, indicating that the scatter of en- 
ergy levels in Fig. [5] is essentially similar to the qh oscil- 
lation across the cut and therefore due to the conformal 
levels spreading in the upward direction. The Coulomb 
qh therefore has the same behavior as the model one, at 



least sufficiently far from the cut. 




FIG. 6. (Color online) Conformal- limit ES of the Coulomb 
v = 1 state for N = 16 bosons and a unit flux added. Delta 
impurity is located at each pole to pin the non-Abelian qh, 
as seen in the e/2 step in the plot of the excess charge (left 
inset). The cut is defined by Na = 8 and I a = 8. Right inset 
shows the zoom on the lowest ES levels. 

We next study the v = 1 Coulomb state, expected to be 
in the universality class of Moore-Read. Abelian vortex 
gives identical counting to the one derived from the root 
partition 0202 ... 02 in Fig. [3] 1 20]. To separate the two 
non-Abelian qhs, one on each pole, we use the method of 
Ref. [H restricting the Hilbert space to the L z = sector 
of the Moore-Read zero-modes. Within this subspace, 
we use a combination of two delta impurities, one on 
each pole, to trap the non-Abelian qhs [io|. To ensure 
that the qhs are indeed separated, we calculate the excess 
charge AQ(0) [2o|. Note that for the fermionic Moore- 
Read state, a two-body potential is necessary in order 
to pin the qhs [l6j- In Fig. [6] we recognize the same 
counting as the one that can be derived from the root 
configuration 1111 ... 11, which confirms that the sector 
change has occurred. 
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FIG. 7. (Color online) Conformal- limit ES of the Laughlin 
model state of N = 12 bosons on the sphere with a unit flux 
removed. The quasielectron is locked to a delta impurity at 
one the poles (m = 0) and displaced half way toward the 
equator (rn — 6). 

Finally, we have also analyzed the ES of the Laughlin 



model quasielectron states, Fig. [7J obtained as ground 
states of the pseudopotential Hamiltonian at one flux re- 
moved compared to the ground state and with a delta 
impurity. Unlike the qh, quasielectron states show a fi- 
nite entanglement gap. This is expected as they are not 
unique and densest zero modes of any local Hamiltonian 
Q . Although the counting of the levels is again a faith- 
ful representation of the edge spectrum, the effect of os- 
cillation across the cut is much more pronounced due to 
their larger size [2l| compared to the corresponding qh 
excitation. This leads to a rapid closing of the entangle- 
ment gap as the quasielectron approaches the equator. 
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